Effective optical constants in stratified metal-dielectric metameterial 
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We present effective optical constants of stratified metal-dielectric metameterial. The effective 
constants are determined by two complex reflectivity method (TCRM) . TCRM reveals full compo- 
nents of effective permittivity and permeability tensors and indicates the remarkable anisotropy of 
metallic and dielectric components below effective plasma frequency. On the other hand, above the 
plasma frequency, one of the effective refractive indexes takes a positive value less than unity and is 
associated with small loss. The photonic states are confirmed by the distribution of electromagnetic 
fields. 



Photonic metamaterials attract great interest as a new 
type of materials including magnetic resonance at optical 
frequency. A famous example exploiting magnetic reso- 
nance is negative refraction^ 

Evaluation of effective optical constants at optical fre- 
quency has been performed for metamaterials of thin 
layers.— In analyzing bulk metamaterials, reflective po- 
larimetry is only way to obtain effective optical constants. 
It is necessary in the analysis to satisfy equation of dis- 
persion. In this Letter we present a way of analysis, two 
complex reflectivity method (TCRM), to determine full 
components of effective tensors of permittivity e and per- 
meability fi. Moreover, the photonic states implied by 
the effective optical constants are examined. 

Figure [T] shows stratified metal-dielectric metamaterial 
(SMDM) and coordinate system. The composite mate- 
rial is obviously uniaxial and the z axis is set parallel 
to the optical axis. We set here metal to be Ag and di- 
electrics MgTV The pair of constituents is selected to 
reduce loss in materials and the interfaces 4 

We assume that effective tensors e and \i are diagonal 
and describe local response in media. Principal axes of 
both tensors are set to be x, y, and z axes in Fig. [IJa). 
In the configuration shown in Figfljb), incident light is 
s polarization (that is, Ei n ||y), and the following relation 
is derived from Maxwell boundary conditions concerning 
bulk material in vacuum (ni = 1, e± = 1, and [i\ = 1): 



k z (9) _ ni cos l-r s {6) 
fi x [ii l + r s (9) 



(1) 



We set E(r, t) cx exp(ik • r — iwt) and k = k/fco where 
ko is wavenumber of light in vacuum. The component k z 
represents the refraction; in particular, k z (0) is refractive 
index along z axis. Complex reflectivity r B was obtained 
by the numerical calculation based on scattering-matrix 
method improved in numerical convergence.^ Optical 
constants of silver were taken from literature,— and the 
refractive index of MgF2 was set 1.38. 

Equation of dispersion under s polarization is 



k z (9) 2 , k x {9f 



Hz 



(2) 



After substituting Eq. ([T]) for Eq. @ , two different an- 
gles enable to evaluate s v /fi x and [i x fi z uniquely. The 
products of e z f[iy and ^ y fi x are obtained by permutat- 
ing configuration (x,y, z) — > (y,z,x). In uniaxial media 
of e x = s y and fi x — /i yi only the two configurations are 
enough to determine all the values of tensors of e and 
\i except for the sign. We call this procedure TCRM. 
This method is generalization of two reflectance method 
(TRM),!£ which is valid for materials of fj, = 1. In TRM, 
conformal mapping between reflectance at two incident 
angles plays a crucial role and retrieves the phase of com- 
plex reflectivity. On the other hand, if it is assumed that 
H 1, we need two complex reflectivity to evaluate e and 
[i tensors. Thus, though TCRM stems from TRM, the 
formalism is quite different. 

It is nontrivial how to determine the sign of /i x . Ac- 
tually we take the sign to satisfy \i x w 1 at off-resonance 




(b) 



(-Hi, e l5 n x ) 




where k x {9) = n\ sin#. 



FIG. 1: (a) Schematic drawing of SMDM and coordinate con- 
figuration. Gray indicates metal (silver) layers and write is 
dielectric (MgF2) layers, (b) Optical configuration of TCRM. 
Incident light is s polarization, that is, Ei n ||j/. 
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FIG. 2: (a) Reflectance spectra. Solid line: under normal 
incidence on xy plane. Dotted line: under nornal incidence 
with Ei n ||z on yz plane, (b) The x component of effective e 
and \x. Upper solid line denotes Re^) and dotted line does 
Im(e x ). Lower solid line shows Re(/Zas) and dashed line does 
Im(fi x ). (c) The z component of effective e and fi. The nota- 
tions are similar to (c): replace x in (c) with z. (d) Refractive 
indexes. Upper solid line depicts Re(n z ) and dotted line does 
Im(riz). Lower solid line represents Ke(n x ) and dashed line 
does Im(n x )- 



energy of 1.5 eV and to connect pb x at resonance with- 
out discontinuous jump. Once the sign of \x x is identi- 
fied, other components of e x , and e z are evaluated 
uniquely; besides, the refractive component k z is also 
evaluated uniquely by Eq. (JT} . As shown in Fig. [21 this 
way determines effective tensors e and /i in a wide energy 
range. 

FigurefS^a) shows reflectance spectra in the two config- 
uration for TCRM; solid line indicates reflectance spectra 
under normal incidence on xy plane and dotted line under 
normal and E; n ||z incidence on yz plane. The periodic- 



ity of SMDM is 75 nm, and thickness of Ag and MgF 2 
are 15 and 60 nm, respectively. Surface layer parallel to 
xy plane is set to be MgF 2 and the thickness is set 30 
nm. Reflectance spectrum shown by the solid line was 
evaluated for thick SMDM of 2000 periods in numeri- 
cal evaluation, and reflectance drawn with dotted line 
was calculated for 50 mm thick SMDM. The thickness 
are enough to ensure that the two reflectance (and r s ) 
are numerically same with infinitely thick SMDM. Nu- 
merical accuracy is tested with changing the number of 
Fourier harmonics employed in the numerical calculation 
for r s and is estimated within a few percent. 

Figures [5Jb) and (c) respectively displays x- and z- 
components of effective e and /i. In Figs.[2]Jb) and (c), up- 
per solid lines denote Re(£j), and lower solid lines stand 
for Re(/ij) (i = x, z). In this TCRM analysis, two angles 
of and 15 degrees are employed. Other pairs of an- 
gles derive same results within the numerical errors. It is 
therefore confirmed that TCRM determines full compo- 
nents of e and \i with satisfying equation of dispersion in 
SMDM. The component e x represents typical permittiv- 
ity of Drude metal below 3.5 eV. The effective plasma fre- 
quency w p>e fT is located at hu> Pje g = 2.6 eV. On the other 
hand, e z implies that SMDM is dielectrics for the light 
of E||z. Below 3.5 eV there exist no prominent magnetic 
resonance in SMDM, while magnetic resonance appears 
at 3.7 eV. The resonance deviates from simple Lorentzian 
dispersion and is associated with resonant behavior of 
permittivity. The 3.7 eV resonance is thus mixed one 
of electromagnetic (EM) components and shows complex 
response due to anisotropy of SMDM. This resonant be- 
havior is discussed later in more details. 

Figure HJd) presents effective refractive indexes. 
SMDM shows a metallic color similar to gold in seeing 
from xy plane, while it is transparent dielectric material 
of refractive index 1.8 in visible range for the light of 
E||z in yz or xz planes. In particular, loss in SMDM is 
quite suppressed in the wide range of 2.6-3.5 eV. This is 
a feature of SMDM. Although the permittivity of bulk 
silver implies large loss in the energy range, the compos- 
ite of silver and MgF 2 is nearly free from loss. This result 
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FIG. 3: (a) Effect of refractive index n z at 3.0 eV. Solid line 
represents the electric field E y (z) of incident light (Ei n ||y, 
normal incidence on xy plane) in vacuum and the refracted 
component in SMDM. Dashed line: the E y (z) evaluated by 
using n z ; Fig. f^d) shows n z = 0.67 + 0.04i at 3.0 eV. 
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FIG. 4: Profile of electric field E y at 3.74 eV (solid line) 
under normal incidence on xy plane. Dashed line is the profile 
calculated from the effective refractive index n z at 3.74 eV 
in Fig. [2jd). Weak reflection is not included in vacuum for 
simplicity. 

suggests that EM fields in metallic layers help transmis- 
sion highly efficient; that is, it suggests that there exists 
EM-field enhancement. Indeed, above w Pj0 ff, several lay- 
ers of metal and dielectrics shows high transmittance. 9 
Figure WLA) indicates that high-efficient transmission in 
SMDM is associated with effective refractive index n z of 
< Rc(n z ) < 1. The effective index thus provides a con- 
cise description for the photonic state. Next, we examine 
EM-field distribution in SMDM to clarify whether the 
state is actually realized. 

Figure [3] displays the E y profile along z axis at 3.0 
eV. Incident light travels along z axis and illuminates 
on yz plane with E; n ||y. Solid line shows incident light 
in vacuum and the refracted component in SMDM. In 
this configuration the effective refractive index is n z = 
0.67 + 0.04i from Fig. f^d). The E y reproduced by n z 
is drawn with dashed line. The wavelength in SMDM is 
obviously close to A/Re(n z ) and longer that A (A: wave- 
length in vacuum). Consequently, it is definitely con- 
firmed that the effective description for SMDM works 
well. In particular, the effective description is excellent 
agreement within a A/{4Re(n z )} scale from the inter- 
face, and suggests that TCRM mainly analyzes the EM 
responses within the A/{4Rc(n z )} depth. The index n z 
represents the low-loss photonic state in a wide range of 
2.6-3.5 eV, where the phase velocity exceeds the velocity 
of light c in vacuum. Figure [3] also shows that the EM 
wave in SMDM is not ideal plane wave and is affected by 



periodic structure. It is one of the inevitable properties 
in mesoscopic metamaterial optics. 

At the end of discussion, we refer to the unusual behav- 
ior of Im(n z ) at 3.7 eV in Fig. HJd). At the resonance, 
Im(n z ) takes negative value and seems to imply expo- 
nential growing wave. However, as shown in Fig. 01 the 
EM distribution at 3.7 eV (solid line) does not show any 
such growing wave, and indicates EM-field enhancement 
in metallic layers. The index n z implies the profile drawn 
with dashed line in Fig. [4J and roughly reproduces the 
EM fields in SMDM within a half wavelength scale from 
the interface. Why does such unusual behavior appear? 
In Fig. [5J the wavelength in vacuum is A = 331 nm and 
the periodicity of SMDM is 75 nm. It therefore seems 
still possible to use the effective description; however, the 
other component is Re^-,,) = 2.0 and A/Re(na;) = 166 
nm which is about twice the periodicity. The condition in 
Fig.[4]is likely close to limits of effective description. Ad- 
ditionally we note that the limits of effective description 
were eagerly debated in GHz range j 1 ^ 11 one of the limits 
is at present understood to be A/|Re(n G ff)| ~ periodicity. 

From numerical results and discussion, we can extract 
a few lessons in TCRM analysis: (i) as for wavelength A 
in vacuum and periodicity a in metamaterial, it is at least 
necessary that A/|Re(n e ff)| > 2a, in order to obtain usual 
effective refractive index, Im(n e ff)> 0. We state that this 
condition is satisfied below 3.5 eV in Fig.[2]and that both 
e and fj, are determined successfully; (ii) effective optical 
constants are determined by the optical response mainly 
within a A/|4Re(n e ff)| scale from the interface. 

In conclusion, we have presented the scheme of TCRM, 
applied TCRM to SMDM, and revealed the full compo- 
nents of e and /i. Analysis of EM distribution has con- 
firmed that the effective refractive index well describes 
the photonic states in SMDM below 3.5 eV. The effec- 
tive index indicates that the low-loss photonic states with 
phase velocity larger than c exists in a wide energy range 
above w PiC ff. 
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